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The dislocation-induced random wnlk of an aluminum-rich droplet migrating in a thermal gradient along 
ihc <100> axis in silicon is shown to lead to an appreciable mean square displacement of the droplet from 
its original <I00> thermal trajectory over migration distances of about 1 cm. When ordered arrays of 
droplets arc misruled to produce columnar diode arrays in silicon, this random walk causes disregistry m 
the arrays, The dtsrr«islry can be minimised by introducing a small ofT-<IOO>-axis thermal gradient 
superimposed on the major <100> thermal gradient. This ofT-axis component of the thermal gradient 
overwhelms dislocation intersection effects which cause random walk and leads to a uniform nonrandom 
displacement which preserves the registry in the columnar diode arrays. 



PACS numbers: 66.30.Kv. 6l.70.Mf, 8l.20.Hy 

t. INTRODUCTION 

The migration of a liquid! droplet in a thermal gradi- 
ent in a solid is caused by the dissolution of atoms of 
the solid on the hot side of the droplet, diffusion of the 
dissolved atoms across the droplet, and the subsequent 
deposition of the dissolved atoms of the solid on the 
cold side of the droplet. The rate of migration may be 
limited eiiiier by the diffusion rate of dissolved solid 
atoms across the droplet or by the atomic dissolution 
and deposition rates at the solid -liquid interfaces of the 
droplet. In this latter case, the migration rate of the 
droplet will be sensitive to the crystal perfection of the 
solid since dislocations and/or grain boundaries which 
intersect the droplet can significantly increase inter- 
face dissolution and deposition rates by providing a 
source of atomic ledges and kinks on the interfaces of 
the droplet. 1 " 3 

A droplet migrating up a thermal gradient will collide 
with dislocations in its path, Such intersecting disloca- 
tions will increase the dissolution or deposition rates 
on those droplet interfaces which they cross. With 
random dislocations, this selective dissolution and 
deposition will occur at random on the various droplet 
facets. For those droplet facets which are not perpen- 
dicular to the thermal gradient, such selective disso- 
lution will cause the droplet to be randomly displaced in 
a direction perpendicular to the thermal gradient. A 
large number of sequential dislocations intersections 
during droplet migration will consequently result in a 
substantial diffusion-type displacement of the droplet 
perpendicular to its original path. 

The random-walk displacement of migrating droplets 
from their thermal gradient trajectories is a technolog- 
ically important problem sir.ee such displacements 
limit the quality of columnar diode arrays, produced by 
migrating a large array of aluminum- rich droplets 
through //-type silicon crystals 4,5 and used in electronic 
devices, radiation detectors and vertical-junction solar 
cells. In this paper, the effects of dislocation density, 
temperature, droplet size, and the crystatlographic 

rs* tsslrrn^i ,K«„„»:«_ ~ - » I. M J - . . --.-II- -1 : - 

v. * u i< iv.k iiiiginwtMi >.. » . t.^imu !.•»• * i<« wyiu - w an\ dis- 
placement of aluminum-rich droplets migrating up a 
thermal gradient in silicon will be examined. It will be 
shown that such displacements can be minimized by the 
selection of crystals with certain dislocation density 



ranges and by the appropriate choice of droplet migra- 
tion directions in the silicon crystal.- 

II. EXPERIMENTAL PROCEDURE 

Previous experience has indicated that the <100) 
direction is the most suitable droplet thermornigration 
direction in silicon because of such factors as droplet 
stability and ease of droplet penetration at the starting 
(100) vacuum-solid interface. 4-0 Consequently, (100) 
l-in.-diam single-crystal ingots of »-type 10-ttcm 
silicon containing between 10 and 10 7 dislocations/rm 2 
were obtained for the experiments. The as-obtained 
crystals were sliced into I -cm-long ingots. The faces 
of these ingots were polished, oxidized, and patterned 
by conventional photolithography and etching techniques 
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FIG. 1. Sequence of schematic diagrams showing the experi- 
ment.-.! prrt.ic.inrc ;:.-.ci to produce an array of' aluminum -rich 
droplets in silicon, (a) Patterning' and etching of an array of 
holes In a silicon matrix by photolithography and etching 
techniques; (b) vapor deposition of ?.n aluminum film In holes; 
(c) removal of excess aluminum by polishing; (d) alloying and 
migration of dioplets in a f ormal gradient. 
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FIG. 2. (a) Wafer face where droplets • niercd. p-type faceted droplet trails In »-typc matrix arc revealed by chemical staining. 
Droplet array has migrated exactly along the (lOO) direction. (lOOx) (b) Wafer face of (a) ingot where droplets exited, Uaculam- 
wulk displacements of droplets perpendicular to the thermal trajectory Ihrmmh the wafer have caused the Initially regular array 
to become irregular. Droplet migration was exactly along n (100) direction and droplet trails arc faceted. (lOOx.) 



to produce a 50x50 square array of holes 30 p deep on 
20 ml centers (Fitf. I). 5 An aluminum film was depos- 
ited into the array of holes. The ingot with the patterned 
array of. aluminum dots was then heated and subjected to 
a thermal gradient of 50°C/cm. The aluminum dots 
alloyed with the silicon and melted to produce an array 
of aluminum-rich droplets which migrated up the ther- 
mal gradient and through the I -cm-thick inpot. *~ n The 
silicon redeposited behind the migrating aluminum-rich 
droplet was doped sufficiently with aluminum to convert 
the originally w-type silicon matrix 'to />-typc silicon in 



the columnar droplet trail. Experiments were carried 
out over a temperature range 945-1330°C. 

The ingot after Ihermomigratlon was ground and pol- 
ished on both the droplet entrance [Fig. 2(a) | and exit 
sides [Fig, 2(b)] and the p-type trails left behind by the 
aluminum-rich droplets in the «-type ingots were 
revealed by conventional chemical staining. Twenty-five 
hundred interncighbor distances, X % between droplet 
trails on the same row for each of the fifty rows were 
measured on both the entrance and exit sides of each 



TABLE I. The mean square displacement along one dimension . dislocation density, temperature, droplet size, droplet morphol- 
ogy, nod droplet migration direction for the droplet arrays examined In the present investigation. Ml droplet arrays have 
migrated 1 cm through the silicon Ingots. 
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FIG. :J. (a) n ml (b) The ontrnncj and exit sides of a wnfer, respectively, where an array of droplets had migrated along but 
sltnliUy o. f f (he (100) nxl n hrcntifte of a sUphtly off-nxis thermal K radlcnt, Droplet trails arc rounded and tear-drop shaped. The 

rwmi#m-wnlk displacements of droplets arc minimised by mipratins slightly off the 000> axis as can be seen by comparing 

with Figs. 2(n) and 2(b), where droplets migrated exactly along the (lOO) direction. (105x.) 



ingot and the mean square displacements (X 2 ) resulting 
from random walk of the droplets away from their 
thermal gradient trajectories were determined from Eq. 
(1). 



, „, U,-U,))» (,v, 

V ' ; ~ 2 2 



(1) 



where A' ( and X f are respectively the interneiplibor 
distances on the initial (entrance) and final (exit) sides 
of the ingot. Been use Die droplet array was deposited on 
the entrance side by photolithography techniques, the 
initial entrance side disregistry of the array { A' ( 
~(A*,>) 2 was negligible and could be ignored in com- 
parison to the exit side disrogtstry (X f -<A',)> 2 re- 
sulting from the random walk of the migrating droplets. 
The denominator of 2 in Eq. (1) accounts for the fact 
that we are measuring. the intemeighbor distance 
between two randomly walking objects and not simply the 
distance of a single randomly walking droplet from a 
iixed reference trajectory. In Table I, the random-walk 
displacement {(X 2 )) 1 ' 2 is listed in order of decreasing 
displacement for the various conditions listed. Also, an 
estimate of the number of displacements ;V experienced 
by a droplet and the average displacement step S was 
obtained by sectioning the sample parallel to the colum- 
nar. droplet trails and examining these doped trails with 
an infrared transmission microscope. In addition to the 
experiments listed in Table I, droplet migration studies 
were also carried out in crystals which had been 
mechanically deformed \% at 1100 °C to produce a high 
dislocation density (~l(J fl dtsl/cm") and in several crys- 
tals with a low dislocation density (10 1 disl/cm 2 ). Dis- 
location densities were measured in the bulk silicon 



crystal after droplet migration with Sirtl etch, a dis- 
location etch for silicon. 

III. EXPERIMENTAL RESULTS 

The most striking experimental result is that droplets 
migrating exactly along the (100) direction suffer rela- 
tively large mean square displacements [Figs. 2(a^ and 
2(b)] as compared to droplets migrating slightly oft the 
(100) axis [Figs. 3(a) and 3(b)]. The doped trails of 
droplets migrating exactly on the (100) are faceted 
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FIG. 4. Longitudinal section of a droplet trail in silicon with a 
low dislocation density of 50 disl/cm 2 . Trails in such material 
arc generally straight. (100x.) 
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FIG. 5. Longitudinal section of n 
droplet trail In fill Icon with a medium 
dislocation density of 10 s dial/cm*. 
Trnlls In such material exhibit many 
relatively large random displacements 
perpendicular to the thermomlgration 
direction. (2Q0x.) 



while those slightly off axis tend to be rounded. 
2 and 3). 



(Figs. 



The effect of dislocations on the mean square dis- 
placement of "ihc droplets varies as the dislocation den- 
sity changes. For low dislocation densities (<5xifj 4 
disl/cm 3 ), the mean square displacements are relati- 
vely small as' might be anticipated. For medium dis- 
location densities (r^lO* to 5'xiO 7 disl/cm 3 ), the mean 
square displacement of the droplets is a maximum. 
Finally, the mean square displacement decreases with 
further increase in dislocation density and becomes 
relatively small again at high dislocation density (10 fl 
disl/cm 2 ). Figures 4. 5, and 6 show respectively, 
metallographic stained longitudinal sections of droplet 
trails in crystals with dislocation densities of «5 x 10 l 
disl/cm 2 , lO* disl/cm 2 , and 10 fl disl/cm 2 corresponding 
to the low, medium, and high dislocation ranges dis- 
cussed above. Both the tow and high dislocation density 
samples have relatively straight trails in comparison 
to the medium dislocation density sample whose trails 
exhibit many random displacements perpendicular to 
the thermal gradient trajectory of the trail. In all cases, 
the dislocation density in the doped silicon trail rede- 
posited by the migrating droplet was found to be much 
less than the dis location density in the surrounding 
silicon matrix. (Figs. 4-6) No dislocations or disloca- 
tion network was found at the doped trail-matrix inter- 
face as shown in Fig. 7, 

No relationship was found between the magnitude of 
the random -walk displacement and the droplet size or 
the Umnorahiro, (Nctc, hrr.vcvcr, that for cquru values 
of(A^>, an optical illusion suggests that arrays of 
smaller droplets are more irregular than arrays of 
larger droplets. ) 



IV. DISCUSSION 

In the discussion we will offer explanations for the 
following observations. 

(1) The large random -walk displacements and faceting 
of droplets migrating along the < 100) direction. 

(2) The smaller random-walk displacements and 
rounding of droplets migrating slightly off the (100) 
direction. 





FIG. 6. Longitudinal section of a droplet trail in silicon with 
n high dislocation density. Trails in such material are 
straight. (2f>0x.) 
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Hon direction will cause a net mean square displace- 
ment of the droplet from its initial theU f j ^ 



Consider a droplet that has already suffered 
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(3) The effects of three ranges of dislocation densities. 

(4) The low dislocation densities in t) 10 doped 

• edepositcd-sitkon trai!s [eti bohind U)e dl . op!ofs 

A. Migration along the (100) direction 

First, wc win bc(tjn by C3lcv}!tUnethc ran(Jom _ 
d isp accments of droplets migrating exactly along the 

tie?iOm e r "h AU,n,inUm - riCh dr °" ,ctS mi e-i ««"»« 
the ( 00) direction in silicon are four-sided pyramids 

wUh e apex of the pyramid pointing forward I (100) 
^ration direction as shown schematically in Fie 
mtlrf? 3C ! Ual , cross Sectio » a'™K a (oh) planc of - some 
mtgra tmg droplets is .shown in Fig. 9. Decause of the 
droplet symmetry, we need only consider random dis- 
placements stoag cither Che (Oh) for ; , 0 < 0I5) (JifWsUun _ 

The random displacements 5 of a droplet perpendicu- 
lar to the 000) dropiet migration results when a droplet 
nutates through a dislocation forest. A dislocation that 
intersects one of the forward faces of ml-ratine ovrl 
midal droplet will cause preferential dissolution of M« 
ace Decause the forward laces of the pyramidal drop- 

S«?h "V* ' t0 th ° <100> direCtion ' Parentis! 
^solution of one forward face will cause a net displace- 
ment of the droplet perpendicular to the <100> migration 

Sfi^.'I!.!-.. 1 ?!- I "« :c :' 0 -. dislocation IntTSons 
17" "l" """""" " ,vv » »*" lynitiutii'iy us « drootet 

migrates through a crystal. The resulting random 
displacements of the droplet perpendicular to its migra- 
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FIG. 10. The mechanism by which a preferential dissolution 
of n forward face induced by nn intersecting dislocation causes 
a ritiiuuiii- tviiik L)'i£*|>liiuc'HiunL 5 in a ui'u|uul pel jiuiaJicuiai lu 
thn Or»w migration direction. 



displacement X (¥ ., because of AT- 1 previous dislocation 
intersections. The net displacement after the Nth 
dislocation intersection is then 



where S is the displacement caused by the iVth disloca- 
tion. Squaring both sides, 

*S= A"J. 1 + 2X / ,-S+S 2 . (2) 

The average random -walk displacement (A'J) caused by . 
A r displacements resulting from <V dislocation intersec- 
tions is found by averaging the terms of Eq. (2). 

<A?v>=<.V}. 1 )+<S a ). (3) 

where the term 2X jV . l "S averages to zero because there 
is no correlation between the net displacement X v ., after 
N - 1 random steps and the Nth displacement S. Iterating 
Eq. (3) (e. g. ; (Xj.,> + (S 2 ); <A'J. 2 ) = 

W. J >+(S»»we find 

(XJ)=MS 2 >, (4) 

where jV is the number of times the droplet has been 
randomly displaced perpendicular to the thermal gra- 
dient by intersecting dislocations and (S 2 ) is the mean 
square average of the variable length displacements S 
induced by the intersecting dislocations. (A 2 ) can thus 
be calculated either from experimental measurements 
of /V and (S*> from metallographic longitudinal sections 
of trails left behind the droplet (see Table II) or by 
calculating /V and (S 2 ) from theory. 

Experimental measurements of N and (S 2 ) result in a 
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C"lci:!r.!crf va!::c r,. f ((A*)) ,/a r.f l,7G v 10" 2 Cm ?.« com- (C^% ta ) ,/2 = 1. 80 xlO" 2 cm 



pared to an observed value of ((A 2 )) 1 ' 2 at the exit sur- 
face of 1.80x10"' cm (Table II) for the first specimen 
in Table I. 
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FIG. U. Two-dimensional projection of pyramidal droplet 
migrating in the (100) direction. In (n) , the dislocation line 
vector Hcs within the apex angle o of the droplet and only 
penetrates one forward face. In (b) , the dislocation line vector 
lies outside the npex angle range o ami thus penetrates two 
opposite forward faces. 

If intersecting dislocations generate the observed 
random -walk displacements of migrating droplets, it. is 
of interest to the technologist to model this phenomenon 
so that the qnnlUalive effects of changes in such para- 
meters as dislocation densities and droplet sizes can be 
predicted. With this limited goal in mind, we will make 
a rough estimation of (A r2 ) from values of /V and (S 2 > 
calculated from a simple two-dimensional model for 
various dislocation densities and droplet sizes. Using 
a two-dimensional model (and thusja_dj:op leLw ith only 
two forward faces), the number of random displace- 
ments N per unit migration length can be found from the 
number of separate times a migrating droplet lias had 
a dislocation intersecting; only one of its forward faces 
during. a migration distance L, The line vector of a 
dislocation which intersects only one forward face of 
the droplet {see Figs. 8 — 10} must be at an angle such 
that it lies in the projected apex angle n of the pyrami- 
dal droplet [Fig. 11(a)). Dislocations lines at other 
angles will Intersect two opposing forward faces and . 
will not lead to a net sideways displacement of the 
droplet. For pyramidal aluminum -rich liquid droplets 
in silicon bound by four (111) planes, the projected apex 
angle is about \u so that only jff/jr - •} of Die dislocations 
that intersect a droplet will intersect only one forward 
face and thus lead to a displacement S of the droplet. 
During a migration distance L, the total number of 
dislocations that intersect a droplet is equal to the area 
2BL swept out along the sides of the migrating droplet 
times the dislocation density p where li is the base widlh 
of the pyramidal droplet. Only about one-third of these 
intersecting dislocations intersect only one forward face 
and thus lead to net displacements in the one dimension 
that we are considering. Hence, the number of displace- 
ment steps N is equal to the density of single-face inter- 
secting dislocations, ip, times the area 2BL swept out by 
along the sides of the migrating droplet for low disloca- 

N^\pBL t pB 2 «l. (5) 



For high dislocation densities pB 3 » 1, dislocations 
will always intersect both opposite forward faces so 
that none of the "**rd faces of the pyramidal droplet 
will dissolve faster than . other. Thus, no displace- 
ments will result so that N will approach aero at high 
dislocation densities, 

# = 0, pB 2 »l. (6) 

Nuw let us calculate the average displacement S 
generated by a dislocation Intersecting a single forward 
face for the case of small droplet sizes (as used In our 
experiments) and low dislocation densities {pB 2 
Figure 12 shows a droplet being displaced sideways by 
a dislocation that intersects only one forward face of 
the droplet. The droplet wilt follow along the dislocation 
line since once the dislocation intersection reaches the 
forward apex of the droplet by sideways migration; 
forward-directed migration other than along the disloca- 
tion line will cause the dislocation intersection to move 
to the face where preferential dissolution on this face 
will cause the apex of the droplet to move back towards 
the dislocation line. Hence the droplet will tend to move 
along the dislocation line until it hits another disloca- 
tion. H If X =p" l/2 is the average spacing between dis- 
locations, then the net displacement S caused by a dis- 
iucaliuit at it* i angle 6 with the thermal gradient is \siny. 
The. average mean square displacement (S 2 ) Is then, 

(S 2 > = (/ [ 8 A 2 s ln*0d0 )( f * 1 0 s in 2 * fi$ )' 1 

= 0.2X 2 = 0.2p-\ (7) 

where we have integrated inly over those dislocations 
Hnes that lay within the apet of the pyramidal droplet 
and thus lead to displacements. 

Substituting in values of -V and (S 2 > from Eqs. (5) and 
(7) into Eq. (4), we find that the mean square displace- 
ment (X 2 ) of a pyramidal liquid droplet with base length 
B migrating up a thermal gradient a distance I through 
a dislocation forest density p is for low dislocation den- 
sities 

(8) 

(X*> = &BL 9 pB 2 «l 

which is independent of dislocation density in this dis- 
location density range. For high dislocation densities, 
we do not have to calculate (S 2 ) since Eqs, (6) and (4) 
already show that 

(,Y*>*0, o£ 2 -I , (9) 

in agreement with our high dislocation density experi- 
ments (see Fig. 6). Unfortunately, most of our experi- 
ments were In a range which was not clearly either a 
very high or very low dislocation density range so that 
one would expect the measured values of QC*) to be in 
between the values given by Eqs. (8) and (9) as is 
observed. 

An important ratio in technology is the ratio of the 
random-walk distance (CO) 1 ' 2 and the penetration dis- 
tance h of the droplet into the material. For example, 
to make high resolution columnar diode arrays in sili- 
con for imaging x rays or infrared, 5 this ratio must be 
small. From Eq. (9), it can be seen that this ratio will 
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FIG. 12. \n) A droplet mi^nttiui; through a diblocntioti forest 
showing the displacement S generated by a dislocation separat- 
ed by a spacing X. <b) A photomicrograph of a displacement 
produced in n droplet trail by dislocations. Note dislocation 
etch pits near the displacement step on either side of the 
droplet trail. (?r,0x.) 



be very small for high dislocation density silicon. For 
relatively low dislocation density silicon on the other- 
hand. Eq. (8) indicates that "Us ratio will be of the 



order {D/lOL) lf7 % Thus to minimize disreglstry in co- 
lumnar diode ^rays in low dislocation silicon, one 
should use the sm«w. ^ droplets possible. 

B. Migration slightly off the (100) axis 

If instead of aligning the thermal gradient parallel to 
the (100) direction, the thermal gradient Is applied in a 
direction a few degrees off the <100) axis, the thermal 
gradient will favor the preferential dissolution of one 
forward face of the droplet because that face will be at 
a higher average temperature than the other forward 
faces. In the silicon-aluminum system as in most other 
systems, the preferential dissolution generated by a 
higher temperature will override the preferential dis- 
solution produced by Intersecting dislocations. Conse- 
quently the off-axis thormal gradient will override and 
prevent the random dislocation-induced sideways dis- 
placements of the droplet perpendicular to the (100) 
direction that occur when the thermal gradient is exactly 
along the (100) direction^ Instead there will bo a steady 
nonrandom displacement of the droplet perpendicular to 
the (100) axis in a direction dictated by the slightly off- 
axis thermal gradient. 

The substitution of a nonrandom displacement for a 
random displacement is technologically fmnnrrant a* U 
allows one to migrate large arrays of droplets through 
silicon without losing the original registry of the array. 
To maximize the preservation of the original registry, 
the off-axis component of the thermal gradient must 
have its sideways component in a (010) direction so that 
two of the four forward faces of the droplet are clearly 
frworcd for dissolution. Figure 13 looking down the (100) 
axis shows a droplet migrating toward the viewer In the 
(100) direction. With the sideways component of the off- 
axis thermal gradient along the (010) direction, faces I 
and n dissolve preferentially and produce a smooth 
nonrandom displacement of the droplet perpendicular to 
the (100) axis. If the off-axis component of thermal gra- 
dient is aligned in the (013) direction (see Fig. 13), 
random displacements would still be induced by dislo- 
cation intersections in the iV directions and would cause 
a disintegration of an initially registered array of mi- 
grating liquid droplets. 

C. Morphology of the droplet trails 

Droplets migrating exactly along the (100) direction 
leave faceted rectangular trail cross sections behind 
them [see Fig. 2(b)]. The reason that these trails are 
generally rectangular and not square is because of the 



FIG. 13. A pyramidal aluml- 
num-rlch droplet migrating 
In a <100) direction In silicon. 
The forward four faces 
formed by (111) planes arc 
shown In a view looking 
down the <100) axis. An off- 
axis component of the ther- 
mal gradient In the <010) 

HlronHrin will «Mn r«*f//o»M- 

walk displacements of the 
droplet perpendicular to the 
(100) axis. 



<010> 
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uneven dissolution of the four forward faces. Ii a dis- 
location threads through two opposing forward faces 
(e.g., faces I and III In Fig. 13), these faces will dis- 
solve faster than the other pair of opposing faces (e.g. , 
faces n and IV In Fl$. 13} and will spread outward, 
causing the base of the droplet to change from a square 
to a rectangle [Fig. 2(b)]. 

The trails of droplets migrating in a slightly off-axis 
(100) direction appear rounded as can be seen in Fig. 
1(b). The off-axis component of the thermal gradient 
was aligned in a (010) direction such that the forward 
two dissolving faces I and II in Fig. 13 remain flat 
because of dissolution faceting. In contrast, the trailing 
faces (faces III and IV in Fig. 13) have become curved 
because of deposition rounding 1 "* 1 and are leaving ho.Mnd 
a curved trail perimeter. The net result is a tear-drop 
trail cross section as can be seen in Fig. 1(b). 

V. DISLOCATION DENSITY IN THE DROPLET TRAIL 

Few if any dislocations are ever formed in the doped 
silicon trails left behind the migrating droplets. Thus 
the dislocation density in the droplet trails is generally 
much less than the dislocation density in the surrounding 
matrix. This result is not altogether unexpected since 
a slow even liquid epitaxy process (the crystal growth 
ccr.Jlwjr.c ir. droplet migration) is known fro hR capable 
of depositing relatively low dislocation density material 
on even a high dislocation density substrate. 
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I4 Thls approximation will overestimate the displacement S 
since the droplet may break away from an Intersecting 
dislocation line before It runs Into another dislocation 
because (1) the dislocation direction may change from a line 
direction that Is effective to a line direction that is Ineffec- 
tive In causing a sideways displacement of the droplet; 

(2) the droplet is in some cases simultaneously pierced by 
two separate opposing dislocations whose effects cancel out; 

(3) the thermal gradient may pull the droplet away from an 
Intersecting dislocation line before any appreciable sideways 
displacement occurs. 
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